Point defects in models of amorphous silicon and their role in structural relaxation 
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We have used tight-binding molecular-dynamics simulations to investigate the role of point de- 
fects (vacancies and interstitials) on structural relaxation in amorphous silicon. Our calculations 
give unambiguous evidence that point defects can be defined in the amorphous solid, showing up 
as anomalies in the valence-charge/Voronoi- volume relation. The changes in the radial distribution 
functions that take place during annealing are shown to be in close agreement with recent, highly- 
accurate x-ray diffraction measurements. Our calculations provide strong evidence that structural 
relaxation in a-Si proceeds by the mutual annihilation of vacancies and interstitials, i.e., local struc- 
tural changes rather than an overall relaxation of the network. 
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The amorphous phase of tetrahedral covalent semicon- 
ductors is usually described iii. terms of an ideal contin- 
uous random network (CRN)Jd which is fully connected 
and perfectly fourfold coordinated; in this model, the 
angular and radial distributions differ only slightly from 
their crystalline counterparts at short range. Yet, the 
atomic structure of these materials — and in particular 
that of the prototypical amorphous silicon (a-Si) — has 
not been completely resolved and the analogy with the 
CRN not fully established. In particular, the presence of 
point defects in the network and their role during struc- 
tural relaxation remains unclear. This is the problem we 
address here. 

A striking similarity has been noted in a-Si between 
the changes induced by annealing (structural relaxation) 
and those associated with thcrjemoval of damage caused 
by radiation in crystalline Si.0 This has led to the sug- 
gestion that point defects such as vacancies and inter- 
stitials (i) do exist in a-Si and (ii) are closely related to 
structural relaxation. Experimental evidence for the ex- 
istence of such defects, and their role in relaxation, i^ 
based on the kinetics ofJieat release upon annealing, 
Mossbauer spectroscopy,^ Cu solubility and diffusion 
and, more receiitly, high-energy X-ray diffraction (XRD) 
measurements.l3 In the latter experiment, the radial dis- 
tribution function (RDF) of pure, ion-implanted a-Si was 
measured with very high precision, both before and af- 
ter annealing. Both samples were found to be under- 
coordinated and equally dense; the average coordination 
number was found to increase from 3.79 to 3.88 upon an- 
nealing. This, as well as other details of the RDFs, was 
interpreted in terms of the mutual annihilation of va- 
cancies and interstitials during relaxation, which brings 
about an increase in the average coordination number 
without affecting the density. 

The concept of point defects provides an attractive per- 
spective for understanding the structure of amorphous 
semiconductors, since it is an extension of well-known 
ideas developed for the corresponding crystalline phases. 



However, some questions remain before this interpreta- 
tion can be fully accepted. In particular, can theory pro- 
vide further evidence for the existence of point defects 
in these materials, and can they be held responsible for 
structural relaxation? What are the structural and elec- 
tronic properties of vacancies and interstitials in amor- 
phous semiconductors? In this Letter, we demonstrate, 
based on quantum-mechanical tight-binding simulations, 
that the interpretation of structural relaxation in a-Si in 
terms of point defects is indeed correct. We propose spe- 
cific criteria for identifying point defects and show that 
the changes in the RDF observed in the XRD experi- 
ments can be reproduced in detail if defects are assumed 
to exist and to annihilate. 

Our study is based on an atomistic description of the 
material. We employed, as a starting point, a CRN model 
containing 216 atoms prepared using a modified Wooten- 
Winer-Weaire algorithm^ and the empirical Stillinger- 
Weber potential. This sample was then relaxed using 
tight-binding molecular dynamicaj (TBMD) at constant 
volume. The TBMD scheme ensures a proper, and "af- 
fordable" , description of the energetics of the material, 
while also providing information on the electronic struc- 
ture; ab initio simulations are out of reach given the size 
of the simulations needed for the present study. The spe- 
cific implementation of TBMD we used is that of Kwon 
et aZjj which is cut-off at 4.18 A, i.e., beyond the second- 
neighbour shell. The TBMD-relaxed sample has a den- 
sity of 0.047 atoms A~'^, is perfectly fourfold coordinated, 
and the width of the angular distribution is a minuscule 
10.30 °. This model will be taken as representing the 
fully-relaxed material. 

Iir order to assess the role of point defects on relax- 
ation, we also prepared "as-made" samples by introduc- 
ing, by hand, a number of defects in the fully-relaxed 
sample. A vacancy is created, simply, by removing one 
atom from its equilibrium position, and an intestitial by 
adding an atom to the CRN. For the latter case, there 
exist many possibilities; we chose configurations closest 
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to the dumbell, as this is the one with the lowest energy 
in the crystalfl After the introduction of each defect, the 
system was fully relaxed. For vacancies, a static relax- 
ation was found to be sufficient. For interstitials, the 
relaxation is sensitive to initial conditions; a simulated 
annealing algorithm was used to scan the energy land- 
scape for a global minimum by quenching slowly from 
300 K to — the temperature was scaled down by a fac- 
tor of 0.997 at every timestep = 1.08 x 10^^^ s. It should 
be stressed that a detailed study of the complete path 
from the as-made to the fully-relaxed material is beyond 
the reach of atomistic simulations. 

Altogether, 40 samples, each containing a single va- 
cancy or a single interstitial (located at random) were 
prepared. In order to demonstrate that defects do exist, 
it is necessary to be able to locate them with reasonable 
confidence in the structure, without a priori knowledge. 
Evidently, no criterion will give a perfect score, as there 
will always be marginal cases due to the continuously 
disordered nature of the material. We have found that 
a criterion based solely on distance — e.g., to identify 
clusters of threefold atoms which may signify the pres- 
ence of a vacancy — are not sufficiently robust for this 
purpose. Nevertheless, it is possible to identify defects 
by looking for correlations between local properties; the 
assumption, which we have verified, is that the presence 
of a defect is felt more strongly at short range. In partic- 
ular, defects are expected to have a sizcjable effect on the 
volume available to neighbouring atomgl3 — the. Voronoi' 
volume — as well as on their valence charges Oll3 

To illustrate this point, we plot, in Figs. 1(a) and (b), 
the Voronoi volumes and valence charges of all atoms for 
the fully-relaxed sample (crosses) as well as two defective 
samples, each containing a single vacancy (squares); as 
will be shown below, these particular samples are typi- 
cal of the set we generated. Modulo some dispersion, a 
clear correlation can be observed, for most atoms, be- 
tween the volume and the valence charge. Both quanti- 
ties are well bounded and the valence charge decreases 
with the Voronoi volume. This can be understood as fol- 
lows: The Voronoi' volume is proportional to the cube of 
the mean nearest-neighbour distance; thus, small Voronoi 
volumes are related to large ion-ion repulsions which can 
only lead to stable structures if screened by large amounts 
of valence electrons. 

For some atoms, however, the correlation fails: some 
points — all of which belonging to the defective sam- 
ples — clearly fall outside the main region. Examination 
of ball-and-stick models of the structures (insets) reveals 
that these atoms are spatially correlated, and in fact sit 
near the sites where the vacancies were created. We note 
that in one case. Fig. 1(b), the "peculiar" atoms do not 
possess a particularly large volumes (two of them are five- 
fold coordinated and the other is fourfold) , while they do 
in the case of Fig. 1(a) (with all four atoms threefold co- 
ordinated); the Voronoi volume, therefore, is not a good 
probe of the presence of vacancies, in contrast to the sit- 
uation in crystalline material. 



Having established that vacancies can be identified in 
a-Si without a priori knowledge of their position, we now 
present a statistical analysis of all 20 samples containing 
a vacancy. Fig. 1(c) shows the correlation between va- 
lence charge, Voronoi volume, and coordination (as de- 
termined by a distance cutoff argument). Again, most 
atoms are found to belong to the main region, but a 
number of data points fall outside of it. Ball-and-stick 
models of the atomic structures reveal, again, that these 
points correspond to atoms surrounding vacancies. The 
particular examples of Fig. 1(a) and 1(b) are thus rep- 
resentative of all samples studied. Similar observations 
apply to interstitial- type defects, which we do not show 
here for lack of space. 

The distribution of points outside the main re- 
gion is found to correlate to coordination: Threefold- 
coordinated atoms (squares) tend to have large volumes 
and exhibit a "bimodal" distribution of valence charge — 
either large or small compared to the atoms in the main 
region. Fivefold atoms (circles), in contrast, occupy rela- 
tively small volumes and are surrounded by low amounts 
of valence charge. The volume and valence charge of 
fourfold-coordinated atoms, finally, present no peculiar- 
ity; their distribution must be taken as inherent to the 
amorphous structure. Thus, while vacancies are pre- 
dominantly associated to threefold-coordinated atoms, 
as is the case in crystalline systems, fivefold-coordinated 
atoms can also signify the presence of a vacancy in the 
amorphous phase. 

We have demonstrated that defects are stable entities 
that can be identified in amorphous structures. We move 
on to assess their role in relaxation by comparing the 
RDFs of models for the as-implanted and the annealed 
states, displayed in Fig. 2, and constructed as follows: We 
first define a local, atom-dependent RDF, Ji{r), which 
is the number of particles per unit length at distance 
r from atom i. This can then be used to compute the 
RDF associated to a particular vacancy, which is done by 
summing the Ji (r) of all atoms located within a distance 
d from the vacancy; this will be labeled Jy{r). Likewise, 
we may define a Jf{r) for interstitials. We may now 
combine Jy and Jf into the total RDF J(r) of a sample 
containing Ny vacancies and A'l interstitials, on average 
2d apart: 



J{r)^7iy4{r)+nijf{r) 



(1) 



where ny = iVv/(iVv + A^i) and ni = Ni/^Ny + Ni). The 
latter expression is formally equivalent to the usual RDF 
for a sample containing a density p of defects with rela- 
tive concentrations of ny vacancies and ni interstitials. 
For proper statistics, Jy and jf were averaged over many 
samples — 20 with a vacancy and 20 with an interstitial. 
There is a single defect per sample and defect-defect in- 
teractions are ignored. Our results are therefore valid in 
the limit of low defect densities, the latter being deter- 
mined by the value of d as p = 1/^nd^. 

We now examine the structural changes brought about 
by the annihilation of defects. We assume here that the 



3 



structure of as-implanted and annealed samples are sim- 
ilar, except for the defect density and relative concen- 
tration. Thus, we can mimic both states of the material 
by simply varying d and ny (or nj). The precise form 
of the evolution from one state to the other can only be 
obtained through a detailed analysis of the atomic mech- 
anisms that take place during annealing, which is beyond 
the reach of current simulations. 

As noted earlier, the XRD measurements of Laaziri et 
alB suggest that the structural relaxation between as- 
implanted and annealed samples proceeds by the annihi- 
lation of point defects, more precisely the diffusion and 
recombination of vacancies and interstitials, which not 
only cause the coordination to increase, but also is such 
that the density is conserved. Thus, it must be the case 
that 

TV^"" = iV^P - 7V™P, (2) 

where A^™^ and N^™^ are the number of vacancies and 
interstitials in the as-implanted sample, respectively, and 
Ny^'^ is the number of vacancies in the annealed mate- 
rial; we assume here that all interstitials annihilate with 
vacancies, iVf™ = 0, which imposes > 0.5. 

Because density, and thus volume, is conserved during 
annealing, it must also be the case that 

where 2(i™P and 2d^'™ are the average distances between 
defects in the as-implanted and annealed sample, respec- 
tively. Combining eqs. ^ and ^ we thus have 

^ann ^ ^imp 3 / -t-iVj 

V A^^P - iVi"°P ^ ' 

The above constraints reduce from 6 to 3 the number 
of parameters needed to calculate the RDF of the as- 
implanted and annealed materials, viz. A'^"p, iVj'™P and 

We now demonstrate that the mutual annihilation of 
vacancies and interstitials are responsible for the changes 
observed during structural relaxation, i.e., increase in 
coordination number of the first and second neighbour 
peaks in the RDF and corresponding decrease of the 
"noise" in-between those peaks. We show in Fig. 2 the 
simulated RDFs for the as-implanted and annealed mod- 
els, using 7i™P = I (and thus n™P = ^) and (i™P = 5.5 
A. The corresponding densities of defects are 3.1 at.% 
for the as- implanted material and 1.0 at.% for the an- 
nealed one, of the same order as those inferred from the 
experiments of Laaziri et alB 

We indeed observe, in Fig. 2, a decrease of the "noise" 
in between the first and second-neighbour peaks and an 
increase in the amplitude of the second peak following the 
annihilation of defects. Quantitatively, the coordination 
of the second peak — computed by fitting the first half to 
a Gaussian — increases by 0.43, in good agreement with 



experiment (0.28) considering the "error bars" inherent 
to both approaches. 

Concerning the first peak, now, an analysis of the as- 
implanted sample shows that the large-r tail is overesti- 
mated when compared to experiment. This is to a large 
extent an artifact of the computational approach, as com- 
plete relaxation of the defects — in particular the intersti- 
tials — is extremely difficult. In order to eliminate these 
artifacts, we follow the same procedure as that used to 
analyze the experimental RDF: we evaluate the coordi- 
nation of the nearest-neighbour shell by fitting the first 
peak to a Gaussian lineshape and integrating. We find 
the coordination to increase by 0.06 from Z = 3.91 in the 
as-implanted sample to Z = 3.97 in the annealed sam- 
ple. This difference is very significant: for Z = 3.91, 
one atom in 11 is, on average, threefold coordinated, 
while the proportion drops to one in 33 for Z = 3.97. 
Taking error bars into account, this change in coordi- 
nation is in excellent agreement with the experimental 
value of 3.88 - 3.79 = 0.09. It should be noted that the 
precise correspondance between experimental and model 
as-implanted samples is impossible to establish; we are 
really only interested in relative changes upon annealing. 

The agreement between model and experiment — and 
thus the interpretation of relaxation in terms of point- 
defect annihilation — is unambiguously established in 
the inset of Fig. 2, where we plot the difference between 
annealed and as-implanted material for both experiment 
(full line) and the present model (dotted line). In or- 
der to set a common reference, the positions of the first 
peak of the RDFs were forced to match. The coinci- 
dence between the two data sets is striking, even though 
the amplitudes of the various peaks is not perfect, due 
to the approximate character of the computational ap- 
proach. The reduction of noise between first- and second- 
neighbour peaks is qualitatively reproduced, and even 
such fine details as the small peaks or shoulders in the 
second-neighbour shell (3.25-4.0 A), indicated by arrows, 
are correctly reproduced. These results provide strong 
evidence that point-defect annihilation is, indeed, the 
mechanism responsible for structural relaxation in amor- 
phous silicon. 

In summary, we have demonstrated, based on TBMD 
calculations of point defects in amorphous silicon, that 
there exists a correlation between the valence charge and 
the Voronoi volume that is broken in the presence of 
defects. Comparison between model as-implanted and 
annealed samples provides strong evidence for the inter- 
pretation of structural relaxation in a-Si in terms of the 
mutual annihilation of vacancies and interstitials, which 
recombine such as to increase the nearest-neighbour co- 
ordination while keeping the density constant. Thus, the 
concept of point defects is relevant to amorphous silicon. 
Annealing evidently proceeds by local structural changes 
rather than an overall relaxation of the network. 
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FIG. 1: Valence charge vs Voronoi volume for samples with 
and without vacancies, (a) and (b) are for two particular 
cases. Crosses are for defect-free samples while squares are 
for samples containing a single vacancy. The ball-and-stick 
models show the spatial correlation which exists between the 
position of the vacancy (black circle) and the positions of the 
atoms exhibiting an anomalous charge-volume relation (grey 
circles). In (c), the corresponding data is collected for all 
20 samples containing a vacancy, and are sorted according to 
coordination, as indicated. The full lines — which are the 
same in all three panels — are a guide to the eye delimiting 
the "normal" region of correlation. 
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FIG. 2: RDFs of as-implaiitcd and annealed models, for a 
relative concentration of vacancies of n™'' = | and d'™^ = 5.5 
A. The inset shows the difference between annealed and as- 
implanted RDFs for both model and experiment. 



